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Abstract 

We consider a class of electromagnetic fields that contains crossed fields combined with longitudinal 
electric and magnetic fields. We study the motion of a classical particle (solutions of the Lorentz 
equations) in such fields. Then, we present an analysis that allows one to decide which fields from 
the class act as a beam guide for charged particles, and we find some time-independent and time- 
dependent configurations with beam guiding properties. We demonstrate that the Klein-Gordon and 
Dirac equations with all the fields from the class can be solved exactly. We study these solutions, 
which were not known before, and prove that they form complete and orthogonal sets of functions. 

1 Introduction 

Relativistic wave equations (Dirac and Klein-Gordon) provide a basis for relativistic quantum mechanics 
and QED of spinor and scalar particles. In relativistic quantum mechanics, solutions of relativistic wave 
equations are referred to as one-particle wave functions of fermions and bosons in external electromagnetic 
fields. In QED, such solutions permit the development of the perturbation expansion known as the Furry 
picture, which incorporates the interaction with the external field exactly, while treating the interaction 
with the quantized electromagnetic field perturbatively [H HI El HI [5]. The most important exact solutions of 
the Klein-Gordon and Dirac equations are: solutions with the Coulomb field, which allow one to construct 
the relativistic theory of atomic spectra [6j, solutions with a uniform magnetic field, which provide the 
basis of synchrotron radiation theory [7], and solutions in the field of a plane wave, which are widely used 
for calculations of quantum effects involving electrons and other elementary particles in laser beams [8]. 
Another physically important class of field configurations (for solving the relativistic wave equations) is 
a superposition of crossed fields and longitudinal fields. Solutions of relativistic equations with fields of 
this type were first studied by Redmond [9]. The Redmond configuration is a plane- wave combined with 
a constant longitudinal magnetic field. The corresponding solutions have wide spread applications, for 
example, in plasma physics [10] and cyclotron resonance In the works |12 [ [T3 t [14] exact solutions of 
the relativistic wave equations with a generalized Redmond configuration (Redmond field plus longitudinal 
electric fields) were found and used to calculate different quantum effects. In [15] the author has presented 
another generalization of a crossed field, which is a particular (the simplest) case of a vortex field [HI fT7] 
(electromagnetic waves with vortices play a central role in singular optics [18J). He studied exact solutions 
of relativistic wave equations in such a field and he has discovered that it can be used to create a beam 
guide for charged particles. 

In the present article we represent and study new solutions of the Klein-Gordon and Dirac equations 
with a new class of fields, which is a combination of crossed and longitudinal electromagnetic fields. For 
the crossed fields E z = H z = 0, E x = H y and E y = —H x , and they depend on the time t and on the 
coordinate z via a light-cone variable £ = ct — z. In the general case, the amplitudes of the crossed fields 
can also contain a linear ^-dependent combination of the coordinates x, y. Thus, we can interpret the 
crossed fields as plane- waves with amplitudes linearly dependent on the coordinates x, y. One ought to say 
that this combination of crossed and longitudinal fields form a class which is described by several arbitrary 
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£ -dependent functions. This combination of electromagnetic fields is physically interesting, since some 
configurations act as beam guides for charged particles in a similar fashion to which the vortex field acts 
in the work [15]. The aforementioned vortex field is a particular case of our beam guiding configurations. 
It is interesting to stress that all other beam-guiding field configurations do not belong to the vortex field 
class. 

The article is organized as follows: first, we describe potentials for the above mentioned combination of 
electromagnetic fields and we study classical particle motion (that is, solutions of the Lorentz equations) in 
such fields. Then, we present an analysis that allows one to decide which fields from the combination act as a 
beam guide for charged particles. We find some time-independent and time-dependent configurations with 
beam guiding properties. Finally, we study solutions of the Klein-Gordon and Dirac equations containing 
all the fields from the combination and we prove that these solutions form complete and orthogonal sets 
of functions. In the Appendix, we place some technical results. 

The electromagnetic fields we consider are defined by the following potential^]: 

A = \ [Q (0 - A] , A 1 = A x = -T x (0 - H(Z)y , 

A 2 = A y = -T 2 (0 + H (0 x, A 3 = A Z = -\ \Q (£) + A] , (1) 

where 

A = R n (0 x 2 + 2R 12 (£) xy + R 22 (£) y 2 , £ - x° - z = ct - z , 

and <?(£), H. (£), T% (£), Rij (0 = Rji(£) , i,j = 1, 2, are arbitrary functions of £. The corresponding 
electromagnetic fields have the form 

E X = H V = T[ (0 +R n (S) x + [R12 (0 + W (0] V, E z = Q' (0 , 

E y = -H x = T' 2 (0 + [R 12 (0 - W (0] x + R 22 (0 y, H z = 2H (£) . (2) 

They consist of crossed fields and longitudinal electric and magnetic fields propagating along the z-axis. 
In the general case, amplitudes of the crossed fields depend linearly on the coordinates x,y. 
The Maxwell current determined by the field Q has the form 

f = ^(p,o,o,P), p = p (0 = Rn (0 + R22 (0 - G" (0 . (3) 
2 Classical motion 

Let us first examine the classical Lorentz equations 

m c 2 u° = e (uE) , m c 2 u = eEw° + e [u x H] , (u ) 2 - u 2 = 1 , (4) 

where 

dx^ 

= — = x^ = (u°,u) , ds 2 = ^ v dx^dx v ', t]^ v = diag (1, -1, -1, -1) . 
And the Hamilton-Jacobi equation is 

(d S+ e -A^ 2 - (vS-- c \f -my = Q , (5) 

where <S is the classical action. From equations (0]) obviously follow the equations for the kinetic momenta 
= m cu^ = (P°,P) : 

m c 2 P° = e (PE) , m c 2 P = eP°E + e [P x H] , (P ) 2 P 2 = m 2 c 2 . (6) 

1 The four-dimensional coordinates of a particle are denoted as x M = (x° = ct, x 1 = x, x 2 = y, x 3 = z), fi = 0,1,2,3, 
where c is the speed of light. Contravariant and covariant four-vectors are often represented in the form 

a M = (a , a 1 ) = (a , a) , a = (a 1 ) , a 1 = a x , a 2 = a y , a 3 = a z , 
£V = rj^a" , a M = rf- v a v . 

Three-vectors are indicated by boldface letters. 
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In particular, from ([6]), with allowance made for ([2]), we easily obtain 

m c 2 P z = e{PE) + e(P°-P z )E z . (7) 
Let us introduce the generalized momenta p^ according to the well-known relations 

Pu = Py. ~ -An, = -<9 M S . (8) 
One can easily prove that the quantity 

A=P°- Pz (9) 

is an integral of motion. Indeed, (JH]) implies 

A = P°~p z = p°-p z + -g(0 ■ (io) 

c 



Hence, we obtain 
Taking into account the obvious relation 



A = P°-P Z + -£?'(?)?' = P°-P Z + -E Z £ . (11) 
c c 



p° - p 

£ = x° - z = u a -u z = - (12) 

TTIqC 

and the equations J6]) and (7]), we find that (fTTjl implies A = 0, which completes the proof. 
Let us introduce the notation 

A = h\, g(0 = -^g(0,p(0 = \-9(0, rn = ^. (13) 

Then lflO|) can be rewritten as 

P°-P z = M?) , (14) 

and 1)121) implies 

m£=p(0=f* = y ^ , (15) 

which relates the proper time and the parameter £. 
In what follows, we denote 

nj (?) - (0 = (0 , *i (0 = ±ft (£) , (t, j = 1, 2) , H (0 = J-H (£) , 5 = ±S . (16) 

Let us also introduce a 2 x 2 symmetric matrix r = r (£) and the two-dimensional columns F = F (£) and 

The complete integral of the Hamilton-Jacobi equations §5§ for the fields <[2j) can be presented as 
S = -i [A (ar° + z) + T] , T = «+/« + x + w + « + X + P + « + v+F + J ( X + X + m 2 ) p- 1 (£) d£ , (18) 
where the 2x2 real symmetric matrix / = / (?), 

f= f hi (?) As (0 

7 ^ /l2 (?) /22 (?) 

and the real two-column x = X (?) satisfy the equations (see Appendix I) 

P (?) [/' (?) + r (?)] - [/ (?) + iH (0 <r 2 ] [/ (0 - iff (0 <r 2 ] = , (19) 
P (?) [X ' (?) + F ' (?)] - [/ (?) + *H (0 <r 2 ] X (?) = . (20) 

Here, 02 is a Pauli matrix. 



Using (fT9|l and l[20|). we can see that the three independent functions (£) provide a solution to a set 
of three first-order non-linear equations, while the two functions \i (£) obey a set of two linear first-order 
inhomogenous equations, where the functions fij (£) are assumed to be known. One should look for a 
particular solution of equations l|19p , and so the general solution of l(20|) has the structure 

X (0 = fciX (1) (0 + k 2 X {2) (0 + X (0 ■ (21) 

Here, \ (£) ls a particular solution for the set of inhomogenous equations (|20l) ; x' 1 ) (£) and \^ (£) provide 
a fundamental system of solutions for the set of homogeneous equations l[20|): fci and k 2 are arbitrary 
constants (two integrals of motion). Thus, the complete integral (fT8| of the Hamilton-Jacobi equations 
© depends on three integrals of motion, A, k\ and k 2 . 

Having at one's disposal solutions of the equations (fl9|) and (|20|l . one can easily find first integrals of 
the Lorentz equations. Using (j4]), with allowance made for lfl5|). one readily obtains a set of equations for 
the coordinates x and y as functions of the variable £ in the following matrix form: 

P (0 v" (0 + V (0 v' (0 - [r (0 + iH' (0 a 2 ] v (£) - (£) v' (£) — F' (£) = . (22) 

One can easily prove that this equation can be integrated once, 

P (0 v' (0 + [/ (0 - iH (0 <7 2 ] v(0+X (0 = (23) 

(see Appendix II). 

Using identity ((6]) for the kinetic momenta, relations (fT4^l and lfT6|) . we find 

(P - P,)(P + P z ) - mgc 2 + P| + P y 2 (P - P Z )(P - P z + 2P Z ) 
= mlc? + Pl + P 2 y =^p 2 (0 [1 + 22 ' (0] = m 2 + p 2 (0 i/+ (0 v> (0 , 

the coordinate z being a function of the variable £. We finally obtain 

2z ' (0 - m 2 p- 2 (£) - v'+ (0 v' (0 + 1 = . (24) 

Expressions ([23]) and lf24|) are first integrals of the Lorentz equations. 



3 Crossed fields and beam guides 

In this section we study a particular case of the field <[2j) in the absence of the longitudinal field, i.e., pure 
crossed-fields, and discuss how these kinds of fields can be used to create a beam guide for charged particles, 
i.e., fields that limit the motion of the charge around some given trajectories. These guides trap the charge 
in a bidimensional plane perpendicular to its trajectory and they are commonly used in many practical 
applications, e.g., quantum computation [19], high resolution spectroscopy [20], non-neutral plasma physics 
|21| . and mass spectroscopy [22]. We will demonstrate that the Lorentz equations for these crossed-fields 
can be reduced to the classical Newton equation with a bidimensional effective potential. As an example we 
discuss a beam guide created by an electromagnetic vortex [15J. Different from the approximated classical 
analogues used to explain the operation of some RF traps [23], the analysis developed here is exact and 
can be used to describe the precise relativistic motion of the charge. 
In the absence of the longitudinal field, we have 

E z = G' (0 = H z = H (0 = . (25) 

As the only influence of the constant Q manifests itself through the z-component of the electric field, we 
can set Q — 0. So our potential ([1]) takes the form 

A = A z = ~A, A x = (0 , A y = -T 2 (£) , 

A = R n (0 x 2 + 2R 12 (0 xy + R 22 (£) y 2 , £ = x° - z . (26) 

Whence, the fields 

E x = H y - T[ (0 + R 1X (0 x + R 12 (£) y , 

E y = —H x = T' 2 (0 + P12 (0 x + R 22 (0 y . (27) 
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For these fields we can identify the integral of motion (|T0|) with the light-front energy E = X/m, and the 
parameter £ is directly proportional to the proper time s, lfl5|) s = m£/A, with A and m given by (fl3)) . 
Substituting the fields f27j) in the Lorentz equation (|4]) we obtain: 

mPi = AP^ = eAEj. = -eX (V±A + d a A ± ) (28) 

where the symbol _L indicates the perpendicular x and y components of the vectors, e.g., = (d x ,d y ). 
We can eliminate the perpendicular components of the potential ([261 using the gauge transformation 

(£, x, y) = (£, x, y) + d M </> (£, x, y) , 

(j)[i,x,y) = xT 1 {C)+y3 : 2{i) ■ 

So, making use of P^_ = fiAx" , the equation lj2"8| becomes 

Ax'_[ = -|Vxi . (29) 

We can identify the above expression with Newton's non-relativistic equation for the two-dimensional 
motion of a particle with effective mass A moving in the effective potential 



(30) 



Therefore, we can find fields that trap a charge in some point of the x, y-plane without explicitly solving 
the Lorentz equations, just by looking for functions A, T\ and T 2 for which the associated potential U is 
capable of limiting the classical motion of a particle of mass A around this point. 

For the special case of a plane-wave, where E = E (£) and H = H (£), we have A = 0, which generates 
the effective potential 

U(Z,x,y) = ^[xF' 1 (0+y^(0] , (31) 

and consequently creates a force HF± — e{T' x ,T'^) that does not depend on the x, y coordinates. So 
although a plane- wave may limit the motion of a charge around some point, it is not possible to fix the 
position of this point only by manipulating the fields. 

3.1 Time-independent fields 

For a time-independent potential, the boundary trajectories can be found by looking for the minima of the 
surface U (x,y). These points can be found using the standard procedure to determine the maxima and 
minima of a function of several variables [24J. A point (xo,yo) will be an extreme if the first derivatives 
dU/dx and dU/dy vanish at this point, and this extreme will be a minimum if the second derivative 
d 2 U/dx 2 and the discriminant D{x,y) are positive at (xo,yo), 

2 d 2 U 



D u y ) = (——\-(—\ >Q — 

' \ dx 2 dy 2 ) \ dxdy J ' dx 



> . (32) 

In the case of a time-independent potential, the expression (|30"j) for U assumes the form 

1 



U(x,y) = ^ 



xC\ +yC 2 -^ (x 2 R n + 2xyR 12 + y 2 R 2 



(33) 



where Cj and Rij = 1,2) are constants. So the condition 11321) implies 



i?n < and RuR 2 2 > (R12) 2 , (34) 



and, consequently, R22 < and det R 7^ 0. Therefore, the minimum for the potential l|33j) . under the above 
restrictions (which is the unique extreme point of U and, consequently, a global minimum), is the point 

_ -R22C1 — -R12C2 _ -R11C2 — -R12C1 
X ° ~ deti? ,Vo ~ deti? ■ 
However, the fields f27j) associated to the potential (|33|) . 

E x =H y = C 1 -x\R n \-y\R 12 \ , 
E y = -H x = C 2 - x \R 12 \ - y \R 22 \ , 

correspond to a problem of a constant charge density in the x, y plane and a constant current density in 
the z direction, which is nonrealistic. 
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3.2 Periodic time-dependent fields 

For a periodic time- dependent potential, and the special case of linear (or quasi-lineaid) systems (which 
includes the case of quadratics potentials i[30|) ). where we have 

V ± i = M (t) x± , 

with M (i) a periodic time-dependent 2x2 matrix, the stability of the potential i|30|) can be studied using 
the Lyapunov criteria [25]. To use these criteria, we first substitute Newton's second order equation l(29|) 
by a pair of first order equations making ui = x^_ , that turns f29|) equivalent to 

n=sv ± ,vx = (£) ,m = ( M \ t) J) ■ (35) 

The motion is called stable around the point (0, 0) if, for every e > we can find a 5 > such that for 
arbitrary initials values Vi (0) with moduli less than S the coordinates Vi (t) remain of moduli less than 
e for all the time t > 0, i.e., the motion is stable if 

Ve > 0, 3S>0:. |V ± (0)| < S |V ± (t)| < e (t > 0) . 

For periodic time-dependent M (t) there always exists a transformation 7?. that leads to a static problem 
E = VJEIRr 1 = const, having the same stability character as S (see [25], Vol. II, p. 119). So, after applying 
this transformation, the stability of the system can be analyzed by finding the roots A& of the characteristic 
polynomial 

det (e - Al) = . (36) 

The system is stable if: 

1. Re(A fc ) < 0, for all X k ; 

2. The pure imaginary characteristic values Re (A&) = (if any such exist) are simple roots. 
If at least one of the above conditions is violated the system will be unstable. 

Let us use the above procedure to analyze the fields generated by the functions lf60|) of the next section, 
for a pure crossed field (H = 0). In this case, we have 

Rn (0 = Ci + C 2 cosw£, i? 22 (£) = Ci - C 2 cosw£ , 
R12 (0 = C2 sinw£ , J 7 ! = T2 = , w, C1.2 = const. 

Substituting these values in ((301 we obtain the effective potential 

tf(£,z,y)=-|r[(x 2 +y 2 )Ci + 0r 2 -y 2 )<72 cos w£ + 2zyC 2 sin w£] . (37) 
Changing to a rotating frame, that is, making the transformation, 

x ± = i?x x , R (0 = f C ° S ^) Sin /J) ) , (38) 
\_ -sm(w4/2) cos(w4/2) J y ' 

the potential l(371) becomes the following time-independent expression: 

U(£ ) y) = -^ ((i 2 + y 2 ) Ci + (i 2 - f) C 2 ) . (39) 

In the rotating frame the equation of motion ([29|l becomes 

i'i = (rMR- 1 -R{R- 1 )"^ x ± -2i?(i?- 1 )'xl , 

M(0 = ^ f ^ + CianfwO \ 

Aft \ C 2 sm(^) Ci - C 2 cos (w£) y v ' 

2 That is, special systems for which the non-linear terms can be neglected, see Chapter XIV of [25J. 



6 



and the expression ([35]) assumes the form 
V', = HV i , V i = 



x 



I 

u± J ' = ~ 1, RMR- 1 - R (R- 1 )" -2R (R- 1 )' 

where, from f38j) and ([4*0]) . we see that 5 is the constant matrix 

I \ = eC, 

ci + c 2 ct 3 + w 2 /4 icr 2 w / ' C% ~ Xh ' 

where at are the Pauli matrices. The roots A& of the characteristic polynomial (|36| are 

1 



Ai = -A 2 = |\/4ci - +4^/ c 2 -w 2 ci , 



A 3 - -A4 - 4 Cl - ^ 2 - 4^/c 2 - c^d . (41) 

Since each eigenvalue has a negative partner, the two Lyapunov conditions will be satisfied only if Re (A&) = 
0. 

The vortex field analyzed in [15] is a special case of (j3T|) for C\ = and C2 = BqU). In this case, the 
constant potential 1(39]) becomes 

U(x,y) = - e -^(^-f) . (42) 

This potential describes the surface of a saddle that rotates (by [38]) in the x, y plane with angular velocity 
lu/2 in time £ (or angular velocity f2/2 = uj\/2m in the proper time s). The classical motion of a particle 
in such a rotating-saddle potential is well known [2S\. However, we were able to obtain some information 
about the trajectories without really solving the equations of motion. 

The condition Re (Afe) = for the eigenvalues (f4*T]l related to the potential ([4*2]) determines the expres- 
sion 

M>^|Bo|. (43) 

This inequality gives us a condition for which the potential lj4"2]) generates bounded trajectories. The above 
result concurs with the condition obtained in [15] by solving Lorentz's equation |(4]) or the one obtained in 
[23] by solving the Newton's equation ([29]) . 

4 Solutions of Klein-Gordon and Dirac equations 

Solutions of the Klein-Gordon equation <I>a, k (£, V: x > y) f° r the fields |(2]), labeled by the three integrals of 
motion A (see [13])) and k = (ki,k 2 ) (see (|2T]0 . read: 



*A, k (£, V,xy)= Nop- 1 / 2 (0 exp (iS) 

A(0=detB(0, 5(0 = 
T] = x° + z = ct + z, (44) 



_/ xi 1} (£) xi 2) (£) 
x? } (0 ^ 2) (0 



where iVo is a normalization factor and (£) (s, s' = 1, 2) are the spinor components (0 introduced 
in (2D), 

This fact can be directly verified by taking into account that the function A (£) obeys the equation 

P-^ = (tr/)A, (46) 
which is a consequence of the uniform set lf20]) . Indeed, the spinors obey the following equation: 

(x (s) ) ' = P" 1 [/ + ^2] x (s) (x (s)+ ) ' = P" 1 X (S)+ [/ - iHa 2 ] . (47) 



7 



The linear independence of the spinors x implies that the matrix B from (|4"4"|l is nonsingular, i.e., 

A = dct B = X « (?) X f (?) - x { 2 ] (?) xl 2) (?) ^ . (48) 
The function A can be easily calculated. One can see that for real spinors x the following relations hold 

A = i X (1)+ <J2X (2) =► A' = l (x (l)+ ) ' °2X {2) + *X (1)+ ^2 (x (2) ) ' • (49) 
Then, using l(4T|) . we find 

A' = [i/a a + H + ia 2 f - H] x {2) = P~\ {1)+ [ifa 3 + ia 2 f] X {2) ■ 

With the help of an evident identity 

if<7 2 + ia 2 f = (tr/)i(T 2 , 

we finally find . 

Solutions of the Dirac equation \E' Aj k (?, f?, a;, y) for the fields in question can be presented in a block 
form by using the two-dimensional Pauli matrices: 

* Al k (?, n, x, y) = Np~ l (?) y^ACl) exp [iS] ( ™ + J ® ^ ) y (0 . (50) 

Here, the two-component spinor y (?) reads 

y (0 = [cosT (0 + ia 3 sinT(?)] , (51) 
where Vb is an arbitrary constant two-component spinor, and we also denote 

T(0 = J ff(0P _1 (0 ^ • ( 52 ) 
The components Fj, i = 1, 2, 3 of the vector f have the form 

Fx = fix (?) x + [/i2 (0 - (0] 2/ + Xi (0 , F 2 = [fu (?) + (0] a: + /m (0 2/ + Xx (?) , F 3 = . 

Therefore, the classical and quantum-mechanical problems are reduced to the solution of the equations 
(PI) and (|20j). 

We will demonstrate that for a complete solution of the problem it is sufficient to find a special particular 
solution of equations (fl9| . 

The non-linear set of equations lfT9|) can be linearized by the following substitution: 



/ (0 = P (?) [cos T (?) + ia 2 sin T (?)] (0 Z' 1 (?) [cos T (?) - 1( t 2 sin T (?)] , (53) 

where Z (?) is a non-degenerate second-order matrix. Using (fl9|) , we find a linear second-order equation 
for the matrix Z (?), 

p 2 (?) (0 + p (0 p' (0 z' (0 + [h 2 (?) - p (0 f (0] ^ (?) = o , 

f (0 = [cos T (0 - ia 2 sin T (?)] r (?) [cos T (?) + ia 2 sin T (?)] . (54) 
A direct calculation yields 

ru (0 = \ru (?) [1 + cos 2T (?)] + ir 22 (0 [1 - cos 2T (?)] - r 12 (?) sin 2T (?) , 
r 12 (0 = fai (0 = n 2 (0 cos 2T (0 + i [m (0 - r 22 (?)] sin 2T (?) , 

r 22 (?) = ^ 22 (0 [1 + cos 2T (0] + \r n (?) [1 - cos 2T (?)] + r 12 (?) sin 2T (?) . (55) 

In order that the matrix / (?) be real and symmetric, one has to look for real solutions of the equations 
(|54l) that obey the subsidiary condition (the symbol ~ stands for transposition) 

J {€) = J (?) > J (?) = z' (?) z 1 (?) ■ (56) 



Such solutions of |H| always exist, since, in accordance with l|5lj) . one easily finds an equation for J (£), 

p 2 (0 J' (0 + p 2 (0 J 2 (0 + p (0 p' (0 ■/ (0 + h 2 (0 - p (0 f (0 = 0. (57) 

Transposing this equation and using the property f (£) = f (£), we find that the equations for J (£) and 
J (£) are the same, which proves the existence of solutions that satisfy f56j) . 

Having at one's disposal a particular solution of equations ([54| that satisfies the condition (|56| . one 
can integrate the equations ff20|) and (|23"]l by quadratures. Indeed, a direct verification shows that the 
expressions 

v = [cosT (0 + *a 2 sinT(0] Z (0 |« - / Z _1 (0 [cosT (0 - *a 2 sinT(£)] X (Op -1 (0 ^ ' 
X- [cosT^+z^sinT^)]^ 1 ^!^-/ Z (£) [cos T (£) - to* sin T (0] F ' (£) dfj , (58) 

obey equations f20|) and (]23]) . respectively, where if is a column with components fci and fe, and «o is a 
constant two-component spinor. 

Motion in the fields of the type (HI) has been studied in previous works. The authors of |26j found 
the symmetry operators for these fields, and the authors of |27[ 128] 129] 130] found solutions for numerous 
specific fields of this kind. However, in all these specific fields there exist certain transformations that 
diagonalize f (£), so that the set of equations (f54| splits into independent linear equations of second order. 

A considerable progress was made in the work [15], which was the first to present exact solutions for 
a specific field that does not admit any transformations leading to the separation of equations l]5lj) into 
independent equations of second order. The author of [H] examined a particular case of the fields ^ with 
the following choice of functions: 



ft (0 =H(0= 9 (0 = 0, r u (0 - ^22 (0 = ccoswf, r 12 (£) = csinc^ 



(59) 



where w and c are some constants. The solutions of the equations in [15] were obtained in a different 
manner from that of the approach of the present work, and they depend essentially on the specific form of 
the fields. 

Let us note that in the present approach there is no necessity to assume Ti (£) = 0, because these 
functions do not enter the set of equations l(5^1) , and so they may be left arbitrary. 

We have succeeded in finding solutions for the fields defined by the following functions: 



5(^=0, H (£) = H = const, m (£) = c\ + C2 cosu;£ , 
^22 (0 = ci — C2 cos w£, ri2 (0 = c 2 sino;£, u, c ij2 = const . 



(60) 



A peculiarity of ([60]) is the presence of a constant and homogenous magnetic field. Expressions ([59]) provide 
a particular case of ([60]) with ff = ci = 0. 

It is easy to see that in the case of the functions {60]) , the equation {54]) can be written in the form 



A 2 Z" + [H 2 - Aci - Ac 2 (crl)] Z = 0, 1 = (sin ft f , 0, cos ft 



(61) 



where T is 



ITT f£\ Hf 



For c 2 = 0, solutions of this equation are known [18] [28, 29, 30J, and therefore we only need to examine 
the case c 2 7^ 0. A direct verification shows that the expressions 



Z n 
Z21 

Z22 



= A 



A 



D 



fix 

a cos cos at 

2 s 

ft x 

a sin — — cos at - 
2 s 

ft x 

(3 sin — - sin /3£ - 



-B 



ft x 

[3 cos — - sin /3£ 



n 
2 

ft 
2 ' 

n 

2 ' 



_ 7 + c 2 

Aft 

7 + c 2 

Aft 

7 + c 2 

Aft 

7 + c 2 

AO 



. fix . 

sin sin af 

2 

ft a; . 
cos — — sin a£ 

cos cos p£ 

• ^ X at 

sin — — cos [it, 



(62) 
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where 



2 

a = 



H 2 - Xci , Q u , 7 2 # a -Aci.fi 



A 2 4 A A 2 

I 2 = c 2 2 + (H 2 - Aci)fi 2 , A, B = const. 



7 

4 A 



give a solution for the equations (|6T|) that satisfies the condition l|56p . The signs of the quantities a, (3, 7 
may be chosen arbitrarily. The quantity 7 is either real or purely imaginary; the quantities a, (5 may be 
complex. For complex a, (3, in view of the linear character of equations (f6T]l . the real and imaginary parts 
of (|62j) separately provide the sought solutions. The expressions |62|) admit a continuous limiting process 
fl — > in case the sign of 7 (7 being real when fi — > 0) is chosen to obey the condition C27 < 0. By 
carrying out this limiting process (and redefining the constants A, B) we find, due to (|62| and fi = 0, that 



z(0 = 



A cos a£ 






B sin /?£ 



of 

a 2 



H 2 - Aci - Ac 2 ) A" 2 
H 2 - A Cl + Ac 2 ) A" 2 



(63) 



The calculation of the quantities / (£) , x (£) ? u (£)> with the help of formulas (|53|) and (|58"1) . is reduced to 
simple algebraic manipulations and integrations of elementary functions, which we omit. 

5 Orthogonality and completeness relations 

In the case under consideration, it is convenient to define inner products for both scalar and spinor wave 
functions on the null-plane £ = const, see for details (HH [32] ■ Such an inner product for scalar wave 
functions is 



($1, $a)e= / (|Q$i(£,»7,a,v)J $ 2(£,r?,ar,j/) 
+$1 (£, »?j ac, J/) 0*2 (£, J/, x, y)} cMxdy , 



(64) 



where 



For spinor wave functions, the inner product on the null-plane has the form 

(« r i,* a )e= / a?, y) *2(-) (C, * s tf) djjdidi/ 

^P ( _ ) * 1 P ( _ ) = i ( l-a 3 ) = i(4 3 

One can verify that scalar wave functions fH|) obey the orthonormality condition 

$a, k) ? = e 5 (A' - A) 5 (fcj - fci) 5 (fc 2 - fc 2 ) , e = p (0 [p (0 |- 1 

— 1/2 

provided iVo = (32-7T 3 ) . First we note that the following relation holds 

Q *a, k (f , V, x,y)=p (0 $a, k (£, ??, x, y) . 
Then, integration over the variable r\ is reduced to the calculation of the integral 



exp 



(A'-A)r; 



dr] = And (A' - A) 



(65) 



(66) 



(67) 



(68) 



Therefore, we can set A' = A in the integral over x, y. The latter is reduced to the following integral 

poo p 00 

J = I dx dy exp 



— 00 j — 00 



g=l, 2 s=l,2 

= An 2 S ( E xl S) (0 (K - h)) S ( E xi s) (0 (*J - fc s ) J • 

\s=l,2 / \a=l, 2 / 



(69) 
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The product of two (^-functions in the right hand side of 1|69|) can be transformed if we take into account 
the following fact: Let a be a nonsingular 2x2 matrix, deta ^ 0, with matrix elements ay . Then the 
relation holds 

S(anzi + a\ 2 z 2 ) 5 (a^izi + a 22 z 2 ) = | deta| _1 (J (zi) S (z 2 ) . (70) 

Setting in l(70|) : ay = Xi (£) ! z± = k[ — ki, z 2 = k 2 — k 2 , we obtain 

J = 4tt 2 |A (0 |-M (fcj - fci) <S (fc 2 - fc 2 ) , (71) 

where A (£) is given by i[4*8|) . Then the result l(66|) follows. 

The constant spinor Vb in solutions <f5l]l is related to an additional (spinning) integral of motion, see 
for details [30, $Tj. So the spinor Vq (and therefore the Dirac wave function as well) depends on a spinning 
quantum number £ = ±1, Vq = Vq ((). It is always possible to choose Vo (£) such that it obeys the 
following relations of orthonormality and completeness: 

v+ (C) v (0 = S c, c ; E y o (0 V (0 = / ■ (72) 

C=±i 

Taking into account (f72| and the relation 

*(-) a, kl c & »7> y) = (32tt 3 )- 1 /2 A Va (£) exp (iS) ( _^ ) V (£) , (73) 

we can verify that the spinor wave functions (|5T|) obey the orthonormality condition 

(*v, k', C *a, k, f > £ = S (V - A) S (k[ - h) 6 (k/ 2 -k 2 )6 cx , (74) 

provided N = (327T 3 )- 1 / 2 . 

The solutions (|44| and (|5Tj) form complete sets of functions on the null-plane £ = const. 
For scalar wave functions (J44J) , we consider the following integral: 

/OO P OO /* OO 

dA / dfci / dfc 2 |p(0l*A,k(€,»/,ar',» , )*A,k(€,»7,a: ) tf) . (75) 

Integrating over the variables ki, k 2 leads us to the integral: 

Mi= / djfei / dfc 2 expfii? (1) fc 1 +ii? (2) fc 2 W4 7 r 2 ( 5('i? (1) ) ( 5('i? (2) ) , 

J —OO J — OO 

i? (s) - xi s) (0 - 1) + x 2 s) (0 (2/ - 2/) , * = 1, 2 . 

This expression has the form lf70|) . where a = S (£) is the transpose of the matrix B in l[44|) . and Z\ — 
x 1 — x, z 2 = y' — y. Thus we obtain: 

M 1 = 4n 2 \A(t)\- 1 5(x' -x)S(y' - y) . (76) 

After that one can easily integrate over A in 1(75]) to get the following completeness relation: 



dA / dA* / dk 2 \p (0 I $ A , k (£, 1/, x', y') $ A) k (£, ry, z, y) 

) J — OO J — OO 

= S (x' - x) 5 (y' ~y)5 (77' - rf) . 

Similar calculations can be performed in the spinor case. Here we have additionally to use the second 
relation ([?2|) to get a completeness relation for the solutions l|5Tj) : 

/oo />oc />oo 

dA / dfci / dk 2 *+ j A k c (£, r/, a;', y') A , k, c K, »?, v) 
. ._ , ! "OO J — DO J — OO 

= S(x' -x)S(y' -y)8tf ~v)P(-) ■ 
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Appendix 

I. Equations (fT9|) and ([201 are obtained as follows. We search for a complete integral of the Hamilton- 
Jacobi equations ((5j) in the form 

S = ~[X(x° + z)+T] , (77) 

with the function T is 

T = hi (0 x 2 + 2/ 12 (£) xy + / 22 (£) y 2 + 2 [ X i (0 + Fx (£)] x + 2 [ X2 (0 + F 2 (£)] y + a (0 . (78) 

Here, (£) , Xi an d a (0 are unknown functions of the variable £. Substituting the expression ([77}. 
with allowance made for ([78} . into equation ([5]), we obtain a quadratic form in x, y, with coefficients being 
functions of £, that must be identically zero, which is only possible when each coefficient is equal to zero. 



Hence, we obtain the following equations: 

P (f'n + ni) - fit - (J it + H) 2 = 1 (79) 

P (/m + ^22) - / 2 2 2 - (/12 - = , (80) 

P (/12 + ri2) - hi (/12 — H) — f 22 (/ 12 + J7) = , (81) 

P ix'i +F[)- fuxi - (/12 + H) X 2 = , (82) 

p(X2+^2)-/22X 2 -(/i2-ff)xi = 0, (83) 

P"'-X?-X2-™ 2 =0 . (84) 



Owing to (f84|) . the expressions (fl8| and ([78} for the function T are identical. It is easy to see that the set 
of equations fT9"ll - l|8Tj) coincides with the matrix equation (fig} , and the set of equations ([82l) - ([83|l has the 
matrix form ([20} . 

II. Let us show that the equation (|22|) is a consequence of the equation (f23|l - To this end, we differentiate 
(j23"ll with respect to the variable ^, 

pv" + p V + [/' - V 2 ] v + [/ - iHa 2 ] v' + X ' = 0. (85) 

However, the left-hand side of ((85} satisfies the identity 

pv" +pV + [/' - iff V 2 ] « + [/ - 2i?cr 2 ] + x' = A + 5 , 
^4 = ?w" + p'u' - [r + iH 'o 2 \ v - 2iHa 2 v' , 

B = [f + r] v+[f + iHa 2 ] v' + X ' ■ (86) 

In the expression for B we now substitute /' from (| 1 Qjl , x' from ([20} , and v' from (|23} , and thus we easily 
find that B = —F' and that A + B is identical to the left-hand side of (|22} . which completes the proof. 
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